Dynamics of spontaneous emission in a single-end photonic waveguide 
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We investigate the spontaneous emission of a two-level system, e.g. an atom or atomlike object, coupled to a 
single-end, i.e., semi-infinite, one-dimensional photonic waveguide such that one end behaves as a perfect mirror 
while light can pass through the opposite end with no back-reflection. Through a quantum microscopic model 
we show that such geometry can cause non-exponential and long-lived atomic decay. Under suitable conditions, 
a bound atom-photon stationary state appears in the atom-mirror interspace so as to trap a considerable amount 
of initial atomic excitation. Yet, this can be released by applying an atomic frequency shift causing a revival of 
photon emission. The resilience of such effects to typical detrimental factors is analyzed. 

PACS numbers: 42.50.-p, 42.50.Nn, 52.25.0s 



I. INTRODUCTION 



A major, if not distinctive, line in quantum electrodynamics 
(QED) is to study how geometric constraints affect the inter- 
action between atomic systems and the electromagnetic (EM) 
field. On the one hand, this can bring a deeper insight into 
the related physics. On the other hand, phenomena that spon- 
taneously do not occur in Nature can become observable this 
way. Spontaneous emission (SE) is an elementary process in 
QED. One normally associates this with an exponential de- 
cay of a quantum emitter (QE) to its ground state accompa- 
nied by an irreversible release of energy to the EM vacuum 
(we will often use the term "atom" to refer to the QE even 
though this needs not be necessarily an actual atom). How- 
ever, free-space SE can be significantly affected - even in 
its qualitative features - by introducing geometric constraints 
forcing the EM field within a certain region of space lU-dt]. 
Cavity-QED has embodied for a long time the traditional test- 
bed for investing such effects. Yet, growing technologic ca- 
pabilities to effectively confine the EM field within less-than- 
three dimensions and make it interact with a small number of 
atoms are opening the door to yet unexplored areas of QED. 
In particular, a variety of experimental implementations of 
one-dimensional (ID) photonic waveguides coupled to few- 
level systems have been developed. These include photonic- 
crystal waveguides with defect cavities |4[], optical or hollow- 
core fibers interacting with atoms JJl], microwave transmission 
lines coupled to superconducting qubits J(J, semiconducting 
(diamond) nanowires with embedded QDs (nitrogen vacan- 
cies) or plasmonic waveguides coupled to QDs or ni- 
trogen vacancies fioll (see Ref. ifTHl for a more comprehen- 
sive list). In most cases, the number (even at the level of a 
single unity) and positioning of such atomlike objects can be 
accurately controlled. Besides major applicative concerns to 
study such ID systems (e.g. some of them can work as highly- 
efficient single-photon sources) these developments are foster- 
ing a renewed interest in their fundamental quantum optical 
properties. Peculiar interference effects, otherwise impossi- 
ble in higher dimensions, can take place such as the ability of 
an atom to perfectly reflect back an impinging resonant pho- 
ton Il2tl which is at the heart of attractive applications such as 
single-photon transistors [[13] and atomic light switches 



To capture certain features of their physical behavior, it of- 
ten suffices to model ID photonic waveguides as endless. In 
reality, of course, one such structure is terminated on both 
sides, each end typically lying at the junction between the 
waveguide itself and a solid-state or air medium. Thereby, 
light impinging on either end always undergoes some par- 
tial back-reflection owing to refractive-index mismatch. Yet, 
mostly prompted by the wish to realize efficient single-photon 
sources, latest technology is now attaining the fabrication of 
quasi semi-infinite, i.e., single-end, ID structures. For in- 
stance, this can be achieved by tapering the waveguide to- 
wards one end so as to make this almost transparent, while 
the opposite end is joined to an opaque medium lalSt]- The 
system thus behaves as being semi-infinite. Equivalently, it 
can be regarded as an infinite waveguide with a perfect mir- 
ror (embodied by the opaque end). Given this state of the art, 
a thorough knowledge of the emission process of an atom in 
such a configuration is topical. While the analogous problem 
in 3D space has been studied extensively @J, first insight into 
the SE of a QE in a semi-infinite ID wave guid e has been ac- 
quired only recently through semiclassical ll 1 6l [l7ll and quan- 
tum models Unlike the 2D or 3D cases, the peculiarity 
of this ID setup is that the entire amount of radiation emitted 
by the atom and back-reflected by the mirror is constrained 
to return to the emitter, hence with a significant chance to re- 
interact with it. As is typical in such circumstances, due to 
multiple reflections, the atom-mirror optical path length then 
becomes crucial and resonances are introduced in the sys- 
tem. This is witnessed by very recent studies (although not 
focusing on SE), where the waveguide termination was shown 
to drastically benefit microwave-single-photon detection lfl9ll . 
atomic inversion schemes (20] and quantum information pro- 
cessing 12 in . 

In Ref. H, through a stationary approach suited to atom- 
photon ID scattering |12] it was shown the emergence of qua- 
sibound states for a single-end waveguide coupled to an atom. 
Indeed, it is known that an atom can behave as a perfect mir- 
ror itself |12lll4| |. hence effective cavities with atomic mirrors 
can be formed l22tl . 

Here, we study the SE of a two-level system coupled to a 
semi-infinite waveguide through a fully quantum model and a 
purely dynamical approach. By deriving the exact time evo- 
lution, we find that the atomic excitation obeys a closed delay 
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FIG. 1 . (Color online) Setup. A semi-infinite waveguide, whose only 
end lies at x = 0, coupled to a two-level system at x = xq. 



differential equation governing the entire SE process. Non- 
exponential emission with long time tales in general occurs. 
This can feature photon-reabsorption signatures and even ex- 
citation trapping, where the latter means that full atomic decay 
to the ground state is inhibited. We rigorously explain such ef- 
fect and illustrate the corresponding output light dynamics. 

The present paper is structured as follows. In SectionHTJ we 
introduce our model and explain the method we used to tackle 
the SE dynamics. Some assumptions that we make are jus- 
tified. In Section |nll after working out the delay differential 
equation governing the atomic excitation time evolution, some 
typical entailed dynamics are shown. In particular, we illus- 
trate the inhibition of a full atomic decay to the ground state. 
In Section [IV] this peculiar effect is demonstrated by work- 
ing out analytically the asymptotic excitation amplitude (i.e., 
its limit at infinite times). We also provide a physical expla- 
nation by rigorously showing that, correspondingly, within the 
mirror-atom interspace a bound state is formed, whose overlap 
with the initial state matches the asymptotic excitation ampli- 
tude. In Section [Vl we illustrate the dynamics of the output 
light exiting the waveguide and show an interesting method to 
induce an emission revival corresponding to a full release of 
the trapped excitation. In Section[VT] we analyze how resilient 
are these phenomena to typical detrimental effects occurring 
in this type of setups. In Section lVTIl we finally draw our con- 
clusions. The work ends with three Appendixes, where some 
technical details are supplied. 



uide terminates at x — only the sine-like modes are to be 
accounted for. Thereby, the atom is dipole-coupled to mode 
k with strength g k oc sin(fcxo). By neglecting counter-rotating 
terms, the Hamiltonian reads 



H — coq \e)(e\+ I dk <y J taja J t+ I dk (g* (T+fljt+H.c.) , 
Jo Jo 



(1) 



where &+ = &_ = \e)(g\ and k c stands for a cut-off wave vector 
depending on the specific waveguide. The total number of 
excitations is conserved since [H, \e){e\ + Jdk a^a^] = 0. As 
we will focus on the atomic SE, the dynamics occurs entirely 
within the one-excitation sector of the Hilbert space. Thus, at 
a time t the wave function is of the form 



m)) = s(t)\e)\V)+\g) dky(k,t)al\0) 



(2) 



where |0) is the field vacuum state, e(f) is the atomic excitation 
probability amplitude and tp(k, t) is the field amplitude in the 
£-space (the normalization condition \e(t)\ 2 + J dk \ip(k, t)\ 2 = 1 
holds). Using Eqs. (Q]) and (|2]) and the bosonic commutation 
rules [a,t, cik'] = d(k-k'), the time-dependent Schrodinger equa- 
tion <9,| V F) = -in I^F) yields the coupled differential equations: 



s(t) = -ia) e(t) -it dk g k <p(k, t) 
Jo 

d,(p(k, t) = -iio k <p(k, t) - ig* k e(t). 



(3) 
(4) 



In line with standard approaches for tackling similar sys- 
tems lfl2ll . we shall make two main assumptions. First, the 
photon dispersion relation can be linearized around the atomic 
frequency as u>o + v(k-ko), where v = daj/d^^ is the 
photon group velocity and ko is such that w ro = ojq. More- 
over, to simplify our calculations we approximate the integral 

rk c poo 

bounds as J Q 1 dk — > J m dk. These approximations are valid 
because we will focus on processes where only a narrow range 
of wave vectors around k-k^ is involved. Hence, wave vec- 
tors which are far from ko (including k < and k > k c that are 
unphysical) have negligible effect. In the following, we will 
set gk = -\ZFv/n sin kxo, where T is the atomic SE rate if the 
waveguide were infinite (no mirror). This assumption will be 
justified a posteriori shortly. 



II. MODEL AND APPROACH 

We consider a ID semi-infinite waveguide along the .Y-axis, 
whose only termination lies at x — 0. The waveguide is cou- 
pled at x — xq to a two-level atom, whose ground and excited 
states \g) and \e) have a frequency splitting ojq. Thus xq is the 
distance between the atom and the waveguide end, the latter 
behaving as a perfect mirror. We sketch the entire setup in 
Fig. 1. The waveguide supports a continuum of electromag- 
netic modes, each with associated wave vector k, frequency 
a>k and annihilation (creation) bosonic operator (a|). In the 
case of an infinite waveguide, for each k > two orthogonal 
standing modes are possible with spatial profiles oc cos(£x) 
and oc sin(fct), respectively. In our case, given that the waveg- 



m. SPONTANEOUS EMISSION DYNAMICS 

Next, we study the system's dynamics when the atom and 
field are initially in \e) and |0), respectively. The initial con- 
ditions thus read e(0) = 1 and ip(k, 0) = for any k. We start 
by removing the central frequency u>o from Eqs. ([3]) and ©, 
via the transformation s(t) -> s(f)e"" J °' ', tp(k, t) — > ip(k, f)e~' w °' '. 
Eq. (01 is thus integrated in terms of the function e(f) as 
<p(k,t) = -i ^Tvjn j^dse iv( - k ~ ko)< - s ~' ) s{s). Replacing this into 
Eq. (fJJ gives 



s(t) 



"Jo 



ds s{s)e 



,-i'uto(j-j) 



dfcsin (kx Q )e 



ivk(s-t) 



(5) 
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FIG. 2. (Color online) Atomic excitation probability P e (t) = |e(r)| 2 
vs. time (in units of 1/F) in the cases Ttj = 2 (a) and Ttj = 0.1 (b) 
and for = oo (i.e., no mirror; solid black line), (f> = 2nn (blue dashed 
line), (p = n/2+2im (red dotted) and (f> = (2n+ l)n (green dash-dotted). 
In (a), only the range t > td is shown: at earlier times the behavior 
does not depend on tp and is the continuation of the black solid line. 



The integral over k is easily calculated as a linear combination 
of 6(s - t±td) and 5{s - t), where td = 2xq/v is the time taken 
by a photon to travel from the atom to the waveguide end and 
back [see Fig. 1]. Once this is used to carry out the integration 
over s, we end up with a delay differential equation for s(t) 
with associated time delay t^. 



s(i) = -$s(i) + $e%(t-t d m-td) 



(6) 



where 0(t) is the Heaviside step function while the phase 
<p = 2koxo is the optical length of twice the atom-mirror path 
[see Fig. 1]. A similar equation was obtained in Ref. 111911 by 
working in real space. The first term on the right-hand side of 
Eq. © describes a standard damping at a rate F. The second 
term, instead, indicates that atomic re-absorption of the emit- 
ted photon can occur at times t > td- At earlier times, such 
re-absorption term is null since the photon has not yet per- 
formed a round trip between the atom and the mirror. Also, it 
vanishes in the limit t d — > oo since the waveguide then effec- 
tively becomes infinite and, as expected, the atom undergoes 
standard, namely fully irreversible, SE at a rate T according 
to |e(f)l 2 = eT Vt . In passing, note that this in fact justifies our 
choice of the factor oc sin(&xo) in gk [see end of previous sec- 
tion]. 

Eq. (|6]l can be solved iteratively by partitioning the time 
axis into intervals of length td- In Fig. 2, we plot the time 
evolution of the atom excitation probability P e (t) = \s(t)\ 2 for 
different values of <p for Ytd = 2 (a) and Ttd = 0.1 (b). In either 
case, the expected purely exponential decay occurring with an 
infinite waveguide (i.e., the no-mirror case) is displayed for 
comparison. Such behavior clearly takes place even in the 
present setup as long as t < td (independently of <p and td). As 
soon as t > td, however, the presence of the mirror starts af- 
fecting the atom in a way that the dynamics is now strongly 
dependent on <p and t d - For Tt d of the order of one, such as 
in Fig. 2(a), the behavior of the atomic population can de- 
viate sensibly from an exponential decay: it exhibits one or 
more peaks of partial atomic re-excitation and, eventually, a 
monotonic decay. The phase <p affects both the positions of 
such re-excitation peaks and the long-time behavior of P e (t) 
[see Fig. 2(a)]. When instead Ytd ^ 1> such as in Fig. 2(b), 
P e (t>td) drops monotonically with the phase <p simply affect- 



ing the atom's average lifetime. Indeed, in such regime the 
solution to Eq. (O can be approximated as (see AppendixlAT) 



e(/)^e 2'0(t d -t) + e 



1+e'* 



2 > 



e{t - t d ) , (7) 



up to an irrelevant phase factor. The corresponding |e(f)l 2 de- 
cays monotonically since |l + e^^|<|l + ^|. 



IV. ATOM-PHOTON BOUND STATE 

An interesting feature emerges for = 0. Fig. 2 indeed 
shows that, regardless of Ttd, such optical path length in- 
hibits a full excitation decay of the atom on the considered 
timescales. Actually, the atom holds a significant amount of 
excitation even in the limit t — » oo. To show this, we take 
the Laplace transform (LT) of Eq. (O and solve the resulting 
algebraic equation. This yields 

r,. 1 Z 7T » (8) 



J<p-st d 



where s(s) is the LT of s(t). According to the final value the- 
orem, e(?^oo) = lim s ^o[se(i)]. We thus find lr23ll 



s(t — >oo) 



= lim[^)]J( 1 + ^r for ^ 2 ^ (9) 
{ for 4>±2nn. 



In the former case note that, in particular, the asymptotic ex- 
citation increases when td is reduced witnessing the crucial 
presence of the mirror. The lower F the more sensitive is the 
increase, i.e., the less uncertain is the atomic emitted light 
wavelength the more pronounced is the effect suggesting an 
interference-like mechanism behind the phenomenon. Such 
inhibition of spontaneous emission can be interpreted as due 
to a destructive interference between the different paths that 
the emitted photon can take to exit the waveguide. This is 
somewhat unusual: one normally expects a two-level atom in 
contact with an electromagnetic reservoir to eventually release 
its excitation altogether. We show next that this may not oc- 
cur here because the atom and traveling photon can form a 
metastable bound state. 

If existent, an atom-photon bound state l^) = e/JeO^ 
Jdlapi,(k) \g) a| |0) must fulfill the normalization condition 
C^bl^b) = 1 and the time-independent Schrodinger equation 
H^i,) = E 1^). The corresponding equations for the atomic 
excitation Sb and the field tfh{k) can be obtained from Eqs. (f3]) 
and (|4|, respectively, through the replacement d, — > -iE as 



E'si, = -\jTv/n J~dk sin(kxo)ipi,(k) 
[E' -v(k-ko)]tpb(k)= ^Tv/tt sin(kxo)sb , 



(10) 
(11) 



where we have set E' = E - a>o- This is solved as 
ipb(k) = -^T/(7iv) sin(kxo)Eb/(k - k ) with k Q = k^+E'/v. 
The squared norm is therefore given by (^bl^b) = [1 + 
T/(nv) Jdksm 2 (kxo)/(k-ko) 2 ]\sb\ 2 . Calculating the integral 



4 



over k through standard contour-integration methods and im- 
posing the normalization condition yields 



|£i| 2 = (l + ^ COS(2£oX ))" 



(12) 



We now replace <fib(k) as given by Eq. (fTTT i in Eq. (TTOb so as to 
end up with E' —-T/2 sin(2&o*o)- It is now easy to show that 
E' = 0, i.e., E = a>o: by taking k — ko in Eq. (fTTT i the left-hand 
side vanishes and hence fulfilling the identity imposes koxo = 
nn (n is an integer). This then brings about k(, = ko,(p- 2nn 
and s(t — > oo) = \sb\ 2 - In conclusion, for <f> = a bound state 
I^Pfc) having energy E = a>o arises. This significantly overlaps 
the initial state |e0), which fully explains why full atomic 
decay to the ground state is inhibited for <p = Inn. As is easily 
checked IT24ll . the amount of atomic excitation corresponding 
to this overlap remains confined within < x < Xo, which 
prevents the atom from undergoing a full decay to the ground 
state. 

This is interpreted as follows. As mentioned, an atom on 
resonance with a traveling photon behaves as a perfect mirror 
I12I , [Till . Thus, just like in a Fabry-Perot interferometer, one 
expects a field standing wave to arise within < x < xq when 
the emitted wavelength matches xq. This agrees with the find- 
ings in Ref. fuHl . which were however derived in the strong 
coupling limit. Significantly, our fc-space approach thus al- 
lows to prove that a bound state is indeed created and work 
out explicitly its exact form. 



V. OUTPUT FIELD DYNAMICS 

So far we have focussed on the atomic excitation dynamics. 
A natural way to experimentally test this is to measure the 
light emitted through the free - i.e., non-reflective - end of the 
waveguide. It is then important to study the entailed dynamics 
of such output light. The real-space field annihilation operator 
at position x> can be expressed as 



C(x) = J- I dk cik sin kx, 



(13) 



where the pre-factor stems from the normalization constraint. 
Once applied to the state in Eq. (O this yields C(x) \*¥(f)) = 
i/f(x, t) \g) |0>, where 



>Kx,t) = 



dk (p(k, f) sin kx 



(14) 



can be interpreted as the real-space field amplitude. The 
square modulus of if/(x, t) can be measured via the statistics 
of the field intensity, which is oc <>P(f)| C\x)C{x) |T(/)> = 
\if/(x, f)| 2 - We assume that a photon detector lies at position 
x — xq + d, where d > is the atom-detector distance. Hence 
(see Appendix iBli. 



(15) 
(16) 



m, t) =yj fee iM ww)- e'Ut' - t d w - m 

= J±e ik ° d [s(t') + iu 0S (t')]6(t'), 




FIG. 3. Output field intensity vs. f = t-d/v in arbitrary units for 
<f> = 2nn under an applied frequency shift A(f) (this is plotted in the 
insets in units of T). (a) Ttj = 2 and a steplike A(r). (b) Ttd = 0.1 and 
a sinusoidal A(f). In either case, a revival of the photon emission oc- 
curs when the frequency shift is switched on. Note that, in line with 
Eq. dl5b . at t = tj the intensity exhibits a discontinuity [particularly 
visible in (a)]. 



where f = t — d/v, and the last equality follows from Eq. ((6]). 
Eq. ( fTSI ) shows that the time evolution of the atomic excita- 
tion can in fact be obtained by integrating the output field am- 
plitude over time (once the time delay d/v is accounted for). 
Clearly, for <p = 2nn, which ensures the formation of the atom- 
photon bound state (see previous section), the atom cannot 
fully decay to the ground state and thus less than one photon 
exits the waveguide on average. An interesting simple method 
exists, though, to force the trapped excitation to be released. 
At a time long enough that the unbound excitation has left the 
waveguide, an atomic frequency shift A is applied (this is rou- 
tinely implemented through local fields). As this changes the 
atomic frequency uq and thereby the corresponding ko, <p is 
modified as well. This suppresses the bound state and neces- 
sarily compels the trapped excitation to leak out as light. In 
Fig. 3, we model the frequency shift switch as a smooth time 
function A(f) [in a way that in Eq. ([TJ a>o — > ojq + A(f), while 
in Eq. © e(f) — > s(t) + /A(f)e(f)] and plot the resulting numer- 
ically computed output field intensity against time. Clearly, 
as soon as A(f) + a spontaneous emission revival takes place 
witnessing that release of the bound state excitation has been 
triggered. For Ttd «: 1 [Fig. 3(b)], the initial emission (when 
A is still zero) is about negligible because in such regime most 
of the energy is trapped within the bound state [cf. Eq. (O and 
Fig. 2(b) for = 0]. Moreover, the system responds quickly to 
applied frequency shifts. Thus in this regime the shape of A(f) 
is closely reflected in the temporal profile of light intensity, as 
shown in Fig. [3b) in a paradigmatic case. Such effect has the 
potential to be harnessed to emit single-photon pulses direc- 
tionally and with controllable temporal profiles. This can be 
of concern to a variety of fields especially in connection with 
quantum information technologies. 



VI. RESILIENCE TO DETRIMENTAL EFFECTS 

To assess the experimental observability of the central phe- 
nomena presented so far, we have refined the model to account 
for detrimental factors. In addition to the waveguide modes, 
we allow for an extra atomic coupling to a reservoir of exter- 
nal non-accessible modes at a rate r ext . Also, we assume that 
the guide is terminated at x — with a non-ideal mirror of re- 
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FIG. 4. Robustness of P e (t) = <|e(0l 2 ) vs. time for (f> = Inn. (a) We set 
(T+T ext )t d = T tot t d = 2, R = 0.95, <5w = r tot and vary the ratio between 
r cxl and T, keeping r tot fixed, (b) We fix Ttj = 0.1, R = l,r cxt = 
and vary 6a>. 



flectivity R < 1 . Moreover, we add phase noise of the atom by 
adding a small white-noise stochastic term to the excited-state 
frequency as a>(t) = coq+Soj 77(f), where 5co is a fluctuation am- 
plitude while rj(t) is a Gaussian-distributed random variable 
such that (77(f)) = and (r](t)r](t')} = 6(t-t') ((■■■} stands for the 
ensemble average). The corresponding excitation probability 
P e (t) = <|e(f)| 2 > in such non-ideal conditions can be predicted 
through a semi-analytical procedure (See AppendixO. 
The first photon-reabsorption peak of P e (f) occurring for 
Ttd ~ 1 [cf. Fig. 2(a)] is rather robust to dissipation into exter- 
nal modes (expected to be the major detrimental factor affect- 
ing such specific feature). As shown in Fig. 4(a), this is still 
visible even with r ext = F (even higher values anyway ensure 
that SE significantly departs from the mirror-less case). As for 
the resilience of the bound-state effects, which are stronger for 
Ttd^- 1 [cf. Fig. 2(b)], Fig. 4(b) shows that in the case of pure 
dephasing a significantly long-lived excitation trapping still 
survives for relatively high 5a>/r ratios (even several units). 



VII. CONCLUSIONS 

We have investigated the time evolution of spontaneous 
emission for a two-level system coupled to a semi-infinite ID 
photonic waveguide. We have shown that the atomic exci- 
tation amplitude obeys an exact delay differential equation. 
According to this, the atomic excitation undergoes a non- 
exponential decay which can exhibit oscillations (a signature 
of partial photon reabsorption) and long time tales. A full 
decay to the ground state is even inhibited when the emitted 
wavelength matches the atom-mirror distance, owing to the 
formation of an atom-photon bound state which we exactly 
derive. The amount of trapped excitation can be substantial. 
Interestingly, it can be released altogether by applying a fre- 
quency shift resulting in a light emission revival. We have as- 
sessed that such phenomena are rather robust against typical 
detrimental effects indicating that an experimental demonstra- 
tion of their key features could be not far-fetched. 
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Appendix A: Derivation of Eq. 10 

For f < td, the delay term in Eq. (O vanishes and thus e(f) = 
e -n/2 p or t > td an( j Yt d <g 1, the time delay td becomes 
the shortest time scale in a way that it can be taken as the 
differential of time. We thus introduce the discrete variable 
and, accordingly, define e„ = s(ntd). Therefore, s (s n+ i - 
Sn)/td, which once replaced in Eq. (|6]l with (e„+i -e„)/frf gives 
the recursion indentity 



1+e ' — 



1 + 



2 > 



(Al) 



Using this along with the matching condition at f = tdB\ 
e(td) — e 2 , we immediately end up with 



Jjd, 

e 2 . 



(A2) 



By combining the functions for f < td and f > td and reintroduc- 
ing the continuous time through n = t/td, we find the Eq. (7) 
of the main text 



r 



e (f) ~ e 2'0(t d -t)+e 2 



t-td 



0(t - t d ) . (A3) 



Appendix B: Derivation of Eq. i fTBl 

We start by recalling that the integration in Eq. (4) of the 
main text for o>* - a>o + v(k - Icq) gives 



,t) = -ig k f dse- ivik - ko)(l - s, s(s), (Bl) 
Jo 



<p(k 



where gk = y/Fu/n sinfcto and the irrelevant phase factor 
e -i&v jj as b een removed from both functions tp(k, f) and e(f). 
Thus, using that the field amplitude in position space is de- 
fined as if/(x, t) = V2/-7T Jdk <p(k, t) sin kx, we find 



x,t) — — / _ f d s fdic s in (kx^) sin (kx)e 
n Jo J 



(B2) 



The integral over k returns a combination of 6 functions, 
which makes particularly straightforward the time integration. 
This yields 

xj,(x, t) = -iy[& [e ik ^- x> £ (t-^)0(x o -x)0(vt-xo+x) 

+e iko(x-x ) s ^_x^o^ x _ Xo)e ^ vt _ x+Xo) 

-e' ioU+A ' o) e(f- ^)6(x+xo)6(vt-x-xo)] . 

(B3) 
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In the special case x — xq + d (d > 0) we have 

^jU) = -»7i e<<W '"* W) [e(f')0(f')-e ,,6 e(f'-f,/)0(r'-f rf )] , 

(B4) 

where t' = t—d/v as in the main text. Eq. (IB4b is equivalent to 
Eq. (13) of the main text, up to an irrelevant phase factor -i. 

Appendix C: Including detrimental effects 

Here, we briefly explain how we have extended our model 
so as to include losses and perform the robustness study in 
Fig. 4 of the main text. For our purposes, it suffices to adopt a 
heuristic reasoning (a more rigorous analysis yields the same 
results). The inclusion of an external reservoir of external 
non-accessible modes is a routine procedure in the literature. 
It simply amounts to adding a term -% i e(f) on the right-hand 
side of Eq. (O, where F ext is the decay rate associated to such 
unwanted modes. The presence of an imperfect mirror with 
R < 1 will instead modify the delay term in Eq. © since the 
atom will re-interact only with the portion of light which is 
reflected. This suggests the substitution e"- 4 — > re"^, where -r 
is the complex probability amplitude for backwards reflection 
from the mirror (|r| 2 = R). For our purposes, the phase of r 
could be absorbed into <p, however a more careful analysis of 
the ID scattering problem yields r = R + i y/R(l - R). Finally, 
as mentioned in the main text, we include extra dephasing of 
the atom by adding a white-noise term to the excited-state fre- 



quency as a>(t) - uq+Su) 77(f), where 5a) quantifies the strength 
of the dephasing and 77(f) is a Gaussian-distributed random 
variable [(77(f)) = and (i](t)?](t')) = <5(t-t')[(---} stands for the 
ensemble average]. In conclusion, Eq. (6) of the main text is 
modified as follows: 

e(t)+iS(OTi(t)sit)=-^^ s{t) + r Y - s{t-t d )6{t-t d ). (CI) 

For completeness, we also mention that a similar analysis can 
be carried out on the output field amplitude, which modifies 
Eq. (13) of the main text as 

m t) = ^ ik '> d [s(t')e(t')-r^e(t'-t d )e(t'-t d )l (C2) 

To obtain each line in Fig. 4 of the main text, we have inte- 
grated Eq. dClt numerically for 100 realizations in the pres- 
ence of simulated white noise, and then averaged over these 
the resulting probabilities P e (t) = |e(f)| 2 . 

Finally, let us stress again that Eqs. dCll ) and iC2i can be ob- 
tained rigorously by modifying the microscopic model given 
by Eq. (1) of the main text. In particular, one has to in- 
clude both sine and cosine waves for each wavevector k, while 
the presence of an imperfect mirror at x = can be mod- 
eled by adding an extra term in the Hamiltonian of the form 
yC f (0)C(0). Here, V = v y/R/(l-R) and C(x) is the field an- 
nihilation operator in real space as introduced in the main text 
[now however owing to the cosine standing modes C(0) ^0]. 
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